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Abstract. One of the greatest problems of primordial inflation is that the inflationary 
space-time is past-incomplete. This is mainly because Einstein's GR suffers from a space-like 
Big Bang singularity. It has recently been shown that ghost-free, non-local higher-derivative 
ultra-violet modifications of Einstein's gravity may be able to resolve the cosmological Big 
Bang singularity via a non-singular bounce. Within the framework of such non-local cosmo- 
logical models, we are going to study both sub- and super-Hubble perturbations around an 
inflationary trajectory which is preceded by the Big Bounce in the past, and demonstrate 
that the inflationary trajectory has an ultra-violet completion and that perturbations do not 
suffer from any pathologies. 
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1 Introduction 

Primordial inflation is one of the most compelling theoretical explanations for creating the 
seed perturbations for the large scale structure of the universe and the temperature anisotropy 
in cosmic microwave background radiation [1]. Since inflation dilutes everything, the end 
of inflation is also responsible for creating the relevant matter required for the Big Bang 
Nucleosynthesis besides matching the perturbations observed in the CMB, see for instance [2, 
3], where inflation has been embedded completely within a supersymmetric Standard Model 
gauge theory which satisfies all the observed criteria, for a recent review on various models 
of inflation, see Ref. [4]. 

In spite of the great successes of inflation, it still requires an ultra-violet (UV) com- 
pletion. As it stands, within Einstein's gravity inflation cannot alleviate the Big Bang 
singularity problem, rather it pushes the singularity backwards in time [5, 6]. In a finite 
time in the past the particle trajectories in an inflationary space-time abruptly ends, and 
therefore inflationary trajectory cannot be made past eternal [7]. Within the framework of 
Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology, this can be addressed by ensur- 
ing the inflationary phase to be preceded either by an "emergent" phase where the space-time 
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asymptotes to a Minkowski space-time in the infinite past [8], or by a non-singular bounce 
where a prior contraction phase gives way to expansion. Unfortunately, within the context 
of General Relativity (GR), as long as the average expansion rate in the past is greater than 
zero, i.e. H av > 0, the standard singularity theorems due to Hawking and Penrose hold, and 
the fluctuations grow as the universe approaches the singularity in the past, which inevitably 
leads to a collapse of the space-time [9] commonly referred to as the Big Bang (see also [7]). 
Thus, in order to provide a geodesic completion to inflation one needs to modify GR. 1 

There have been some recent progress in ameliorating the UV properties of gravity 
in the context of a class of "mildly" non-local 2 higher derivative theory of gravity [12- 
14], involving terms that are quadratic in the Reimann curvature but contains all orders 
in higher derivatives. It was found that such an extension of Einstein's gravity renders 
gravity ghost-free and asymptotically-free in the 4 dimensional Minkowski space-time. 3 As a 
result the ordinary massless graviton sees a modified interaction and propagation at higher 
energies, but at low energies in the infra-red (IR) it exactly yields the properties of Einstein's 
gravity. It has been demonstrated that such a modification of Einstein's gravity ameliorates 
the UV behaviour of the Newtonian potential, and also permits symmetric non-singular 
oscillating solutions around the Minkowski space-time [14] signalling that the action can 
resolve the cosmological singularities, see also [12, 13, 19]. Similar non-local higher derivative 
theories of gravity have also been considered in connection with black holes [20] , inflationary 
cosmology [21], string-gas cosmology [22] and in efforts to understand the quantum nature 
of gravity [23]. 4 There was also a proposal to solve the cosmological constant problem by a 
non-local modification of gravity [26] (see also [27]). 

What is rather intriguing is that non-local actions with an infinite series of higher 
derivative terms as a series expansion in a' , the string tension, have also appeared in stringy 
literature. One mostly finds them in the open string sector, notably in toy models such as 
p-adic strings [28], zeta-strings [29], strings on random lattice [30, 31], and in the open string 
field theory action for the tachyon [32]. For more general aspects of string field theory see 
for instance [33]. There have also been progress in understanding quantum properties of 
such theories, i.e. how to obtain loop expansion in the string coupling constant g s , which has 
lead to some surprising insights into stringy phenomena such as thermal duality [34], phase 
transitions [35] and Regge behavior [31]. It is indeed a wishful thinking to derive such an 
action in the closed string sector involving gravity. 5 

In Ref. [12], with the means of an exact solution, it was shown how a subset of the 
actions that have been proposed more recently in Ref. [14], can ameliorate the Big Bang 
cosmological singularity of the FLRW metric. The action discussed in [12] only contained 

For general reviews of modified gravity cosmological models see for example [10]. 

2 By non-local we only mean that the Lagrangian contains an infinite series of derivatives and hence can 
only be determined if one knows the function around a finite region surrounding the space-time point in 
question. This property also manifests in the way one is often able to replace the infinite differential equation 
of motion with an integral equation [11]. The theories however are only mildly non-local in the sense that the 
physical variables and observables are local, unlike some other theories of quantum gravity /field theory whose 
fundamental variables are themselves non-local, being defined via integrals. 

3 For attempts along these lines with a finite set of higher derivative terms, see [15-17]. However in four 
dimensions they were not successful in avoiding ghosts and obtaining asymptotic freedom simultaneously, due 
to essentially the Ostrogradski constraint [18]. 

4 For investigations along these lines in the more broader class of non-local gravity theories which include 
non-analytic operators such as 1/D, see [24, 25] and references therein. 

5 We note that significant progress has been made in studying non-local stringy scalar field models coupled 
minimally and non-minimally to the Einstein gravity [36]-[39]. 
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the Ricci scalar and it's derivatives up to arbitrary orders and yielded non-singular bouncing 
solution of hyperbolic cosine type. However, it was not clear how generic these solutions were, 
whether these solutions were stable under small perturbations, and whether these theories 
contain other singular solutions or not. 

Some of these issues, in particular concerning only time dependent fluctuations around a 
bouncing solution were addressed in Ref. [13], where the authors have analysed the stability 
of the background trajectory for very long wavelength (super-Hubble) perturbations, either 
when the cosmological constant is positive, i.e. A > 0, or negative, i.e. A < 0. 6 The latter 
provided a way to construct cyclic cosmologies, and in particular cyclic inflationary scenarios 
that have been studied in [40], and in [41] where cosmological imprints were also discussed. 

In this paper we aim to focus on the case when A > 0. This particular case is very 
interesting from the point of view of inflation, as it paves the way for a geodesically complete 
paradigm of inflation in past and future. It was found that there exists a de Sitter solution 
as an attractor for large enough times before and after the bounce. The weakening of gravity 
in the UV regime facilitates the bounce for the FLRW metric. Our main goal will be to 
study what happens to the non-singular background solution under arbitrary classical and 
quantum fluctuations (super and sub-Hubble type) during the bounce period. We have a 
two- fold aim. First, we want to verify that at late times there are no fluctuations that are 
growing during inflation. This would clearly signal an instability in the system which would 
certainly jeopardize the inflationary mechanism of creating superHubble fluctuations which 
becomes a constant at late times. Second, we want to check that the bounce mechanism itself 
is robust, i.e. the perturbations do not grow unboundedly near the bounce. In the process of 
our investigations, we will develop techniques to deal with non-local perturbation equations 
that one encounters in these theories, and this we believe is a major achievement of the paper. 
We hope that this will help us in analyzing other cosmological issues in future. Finally, let 
us comment on a rather interesting feature we noted. It turns out that the requirement of 
the bounce automatically introduces an extra scalar degree of freedom, and if the mass of 
this scalar is light, it may be able to play the role of a curvaton [4]. We however leave a more 
detailed phenomenological investigation of such a scenario for future. 

The paper is organized as follows: In section 2, we review some of the results for non-local 
gravity models and generalize some of the techniques previously used to find exact solutions. 
In section 3, we derive the non-local perturbation equations for the Bardeen potentials around 
a large class of cosmological background solutions. Next, in section 4 we proceed to solve these 
equations around the "inflationary" de Sitter late-time attractor solutions that exist in these 
models. In particular, we obtain the criteria for which the perturbations freeze out or decay 
after crossing the Hubble radius which allows one to maintain the inflationary mechanism 
of producing scale-invariant CMB fluctuations. Thereafter, we look at the stability of the 
fluctuations in section 5 around the bounce. Finally in section 6, we conclude by summarizing 
our results and discussing their implications for inflationary cosmology. 

2 String-inspired Non-local Gravity 

In this section we begin by introducing the String theory motivated non-local gravity models. 
We will then generalize the methods previously employed in [12, 13] to obtain exact solutions 
in these theories. We will explore the various non-singular cosmologies that are present 

6 Our preliminary investigations indicate that when A = 0, the higher derivative modifications allows one 
to realize the emergent universe scenario. 
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in these theories and in particular focus on solutions which can geodesically complete the 
inflationary space-time. Although the inflaton potential energy does not remain a strict 
constant and changes slowly, for technical reasons in our analysis we will treat the potential 
energy as a "cosmological" constant. This should not however affect the issues we plan to 
address in this paper. 

2.1 Action and Equations of Motion 

In Refs. [12, 13] the cosmology of a special class of ghost free actions was investigated. The 
action has the form 

S = J d 4 xV=g {^R + ^RF{U/M*)R - A + £ M V (2.1) 

where Mp is the Planck mass: Mp = 1/(8ttGn), Gn is the Newtonian gravitational constant, 
A is the cosmological constant, M* is the mass scale at which the higher derivative terms in 
the action become important, and A is essentially a book keeping device to help us keep track 
of the higher derivative corrections. Cm is the matter Lagrangian. We use the convention 
where the metric g has the signature (—,+,+,+). 

The analytic function F(U/M%) = £ f n U n is an ingredient inspired by String The- 

ory as an expansion in the string tension (a'). We recall that the action of the covariant 
d'Alembertian on a scalar is given by 

V^g 

where V M is the covariant derivative. Variation of action Eq. (2.1) yields the following 
equations: 



, oo n— 1 

A 



[Mp + 2\F(n)R]GZ = - ^ fn Yj [g^dp^Rd^-^R + g^d^RdpU^^R 



n=l 1=0 



- 5£ (g^dp^Rd^-^R + ^RU^R)^ - ^RT(n)R5% 



(2.2) 



where Tm^ is the energy-momentum tensor of matter and 

Gt = R>t- \^R , (2.3) 

is the Einstein tensor. We assume that the matter stress energy tensor obeys the usual con- 
servation equation, V^Tm^ = 0, to ensure consistency of the modified Einstein's equations, 
which we will refer to in future as Einstein-Schmidt (ES) equations 7 . The ES equations can 
be written in a reasonably compact form: 

[Mp + 2\T(n)R]GZ = T M £ - A^ + 

+A/C£ - ±8£ (K% + Kx) - ^RF(n)R5£ + 2X(g^V p d u - 8^U)F(U)R , ^ 



7 The generalization of the Einstein's equations to include terms such as EL\ n R was analyzed by 
Schmidt [42], having an infinite series of such terms give rise to Eq. (2.1). 
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where we have introduced two additional quantities: 

oo n— 1 



oo n—1 



K = 9^ /» E ^□'^□ n ~'~ 1 fl , £1 = ^ /« E O l BD n - l R . (2.5) 



n=l Z=0 



n=l Z=0 



Making any progress with the ES equations may seem daunting at first, but previous studies 
have shown that the trace equation can be tractable since all we have to deal with is the 
Ricci scalar: 



oo n—1 



M 2 P R- Xj2fnY,l d ^ nlRdflnn ^ lR + 2nlRnn ' lR ) - 6AQF(D)i? = 4A - T M £ , (2.6) 



(2.7) 



n=l 1=0 

or, in the more compact notation 

MpR - A/C£ - 2\Ki - %XUF{U)R = 4A - T M £ , 

2.2 General Ansatz for finding Exact Solutions 

It has been shown in Refs. [12, 13, 43] that the following ansatz 

DR -nR-r 2 = and n ^ , (2.8) 
is useful in finding exact solutions. From Eq. (2.8) one gets the following relationships: 



U n R = r\ l [R + 



V2 

n 



for n > 0, 



r 2 , 



(2.9) 



F{U)R = F X R + —{T x -f ) where Ji=^(n). 
If the scalar curvature R satisfies (2.8), then equations (2.2) have the following form: 



Mp + 2A ( T(n)R + -(Hn) - ti 



ri 



d^Rd u R 



g ap d a Rd p R + n ( R + 



(2.10) 



+2AJ\ [D»d v R - 6%( ri R + r 2 )\ - A 



52 



F\R 2 §(Ji — fo) 



r 



where J-^ is the first derivative of J- with respect to the argument, and the trace equation 
(2.6) becomes especially simple: 



where 



A X R - XT m (rt) {2 n R 2 + d^Rd^R) + A 2 = -T M £ , 



A\ = Mp — A U«F«(ri)r 2 - 2^(Ji - /o) + 6Jir^ and 



(2.11) 



A 2 = - 4A - A- ( 2J-( 1 )(n)r 2 - 2— (Ji - / ) + 6Jin 



»"2 
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We proceed to consider a traceless radiation along with a cosmological constant. This will 
simplify many of our calculations. Then the simplest way to obtain a solution to Eq. (2.11) 
with Tm^ = is to put A\ = A2 = 0, and impose 



jpW(n) = . 



(2.12) 



This implies 



r 2 = - 



n[M 2 P - QXFm] 
2A[Ji - M 



A 



r 2 M 2 p 



4n 



(2.13) 



Thus we see that within the context of the ansatz, the ES equations and the corresponding 
solutions are characterized by three independent parameters, r%, A/o, and AJ 7 !. Although A 
is only a book-keeping device, and one can set it to one if desired, it is interesting to look at 
the analytic properties of the different parameters with respect to A since the A — > limit 
can mimic the large time limit when one expects to recover GR and all the higher derivatives 
become small. Now, r\ is determined by (2.12) and is therefore independent of A, but ri 
depends non-analytically on A. 8 This could seem to be a potential problem in recovering 
General Relativistic solutions. As we will see later in explicit examples, although ri may 
have non-analytic dependence, it is really just an intermediate variable of convenience, and 
the solutions for the scale- factor need not have any singular dependence in A (see section 2.4). 
This is why one is able to recover the GR solutions in appropriate limits. 

It is important to note that by virtue of the Bianchi identities, solving the trace of the 
ES equations essentially guarantee that we have a solution for the complete ES equations, 
which now simplifies to 



In general one is required to include radiative sources to get exact solutions [13], and we 
need to make sure that the radiative energy density if present is positive and not ghost-like. 
We will discuss this point in the specific context of cosmological solutions in the following 
subsection. 

Let us emphasize that equations up to this point are general and do not take into account 
the properties of the metric. Thus the above procedure can be used to obtain metrics which 
are cosmological (time dependent solutions) in nature or which have spatial dependence. This 
is an important generalization on the analysis discussed in [13], and a major step forward in 
investigating these non-local models. 

2.3 Cosmological Solutions, Bounce and Stability 

We are now going to focus on spatially flat FLRW type cosmological solutions that was 
obtained in [12, 13] in the presence of a cosmological constant and radiation energy density, 



8 Actually, not only T2, A also acquires a non-analytic dependence through the second equation in (2.13). 
In fact, this has to be the case, otherwise we could take the A — ¥ limit smoothly and conclude that GR has 
non-singular bouncing solutions! 






(2.15) 
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with a traceless energy momentum tensor. In the above expression a is the scale factor and 
ao refers to its value at the bounce. In the cosmological context the ansatz Eq. (2.8) becomes 
a third order differential equation for the Hubble parameter, H = a/a, 



H + 7HH + AH 2 + 12H 2 H 



2 ri H 



2 



[2 

6 



(2.16) 



where dot defines derivative w.r.t physical time. The solutions to the above equation, their 
stability, and asymptotic behaviours have been studied in details in [13]. It was found that 
several non-singular solutions exist where the universe bounces from a phase of contraction 
to a phase of expansion, thus providing a resolution to the Big Bang singularity problem. By 
looking at the '00' component of the ES equations, one can calculate the radiation energy 
density at the bounce point, po, which must be positive: 



plays the role of an 'initial condition'. There are three cosmological scenarios one can envisage 
from Eq. (2.16) : 

• A < 0, n > T2 > 0: One obtains a cyclic universe scenario where the universe 
undergoes successive bounces and turnarounds. Qualitatively speaking the bounce is 
caused by the higher derivative modifications, while the turnarounds occur when the 
radiation energy density cancels the negative cosmological constant. Such scenarios 
can give rise to "cyclic inflationary" models if one is able to incorporate entropy pro- 
duction [40, 41], but we will not consider these scenarios here any further. 

• A>0, ri<0=^r2>0: One obtains a bouncing universe scenario where a phase of 
contraction gives way to a phase of super-inflating expansion phase. An exact example 
is the Ruzmaikin-type solution [44, 45], which incidentally contains no radiation [43]. 
Unfortunately, in these scenarios the universe is stuck in a quantum gravitational phase 
accelerating faster and faster, and we never recover a GR phase. Clearly such solutions 
are not phenomenologically viable and hence we disregard these scenarios. 

• A>0, ri>0=^r2<0: This represents a geodesic completion of an inflationary space- 
time via a non-singular bounce. It is easy to see that the action Eq. (2.1) without matter 
admits constant curvature vacuum solutions, the de Sitter solution or the Minkowski 
solution, depending upon the value of the cosmological constant. These background 
solutions are in fact exactly the same as that of the general relativistic de Sitter phase, 
as for constant R, all the higher derivative terms in the field equations vanish, and one 
just obtains from Eq. (2.6) 



This is precisely the de Sitter solution the space-time asymptotes to in the more gen- 
eral case, as the radiation dilutes away very quickly during the cosmological constant 
dominated phase. 

It is this last case that is relevant for inflationary cosmology, and is the main focus of 
our paper. The higher derivative extension of gravity provides us with a plausible answer to 




RdS 



M 2 ' 
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the question, "what happened before inflation?". The general picture that emerges is that 
in the presence of a positive cosmological constant, if a small causal patch is approximately 
homogeneous and isotropic at Planckian densities, the universe would undergo a non-singular 
bounce, following which the universe will start to inflate quickly diluting any matter (radi- 
ation) present in the system, and we essentially will have a pure de Sitter. Of course, to 
be consistent with the inflationary phenomenology, A cannot be a strict constant but has to 
evolve, possibly with a scalar field, the inflaton [4]. For technical reasons we haven't been 
able to incorporate this dynamics yet, but our main aim here is to assess the robustness of 
the bounce mechanism and the inflationary scenario in terms of perturbations. 

What we expect is that with the dilution of radiation energy density during the vacuum 
energy dominated phase, any fluctuations present in the fluid would also get diluted away, so 
that the fluctuations can become dominated by the quantum fluctuations of the inflaton field 
which we know can produce the observed approximately scale- invariant spectrum in CMB. 
Since we do not have any inflaton in our scenario, what we want to check is whether there 
are any growing or unstable modes of fluctuations around the geodesically complete de Sitter 
space-time. If present, this would signal either the fallibility of the bounce mechanism, or 
the breakdown of the inflationary mechanism of generating fluctuations. 

2.4 The Hyperbolic Cosine Bounce 

It is rather fortunate that one can actually find an exact solution for the higher derivative 
extension of gravity under consideration, which serves as a good illustration of a geodesic 
complete space-time. It is given by 

a(t) = a cosh yj^-t, => H = y^ tanh (^2 *) ( 2 ' 19 ) 
This solution satisfies the ansatz with the specific parameter combination 

UR = riR-6rf. 
Substituting T2 = —§r\ into Eqs. (2.13) and (2.17) we get 

A = ^nM| and p = -yAJirf. (2.20) 

Note that the radiation is non-ghost like, provided 

T(n) < 0. (2.21) 

This turns out to be a general constraint which must be satisfied in order to realize geodesi- 
cally complete de Sitter space-time, and we are going to assume that this constraint is sat- 
isfied. The above specific example (2.19) serves as a nice illustration of some of our general 
results and we will fall back on it as appropriate in the remainder of the paper. 

Finally, one may wonder about the special case when T\ = 0, which means there is no 
radiation. Although this is under current investigations, indications are that such theories 
will suffer from the presence of ghosts. Effectively there is a double pole in the propagator 
at p 2 = r\. 



-8- 



3 Perturbations 



In the previous section we have provided a general algorithm for finding cosmological back- 
ground solutions in non-local gravity models which in particular can realize inflationary space 
time at late times. We have also provided concrete examples of such solutions. What we 
now want to address is the robustness of the paradigm by studying the fluctuations around 
the background. This will let us test whether the bounce mechanism is stable, as well as 
whether the inflationary mechanism of generating scale-invariant perturbations can survive 
in the presence of higher derivative modifications. In this section we will provide a general 
framework for dealing with the scalar modes and obtain the coupled (non-local) differential 
equations describing the evolution of the Bardeen potentials. Since this is a somewhat tech- 
nical section, we have provided a summary subsection 3.4, where we enumerate the results 
in a self-contained manner. 

3.1 Perturbations on higher derivative action of gravity 

To compute perturbation equations we use the fact that the background configuration obeys 
the ansatz (2.8) and conditions (2.13). Let us introduce the following notations to separate 
the background and the perturbations: 

R = R B + SR, n = D B + SD, (3.1) 

where the subscript B stands for background. Also, let us denote the variation of the ansatz 
(2.8) as follows: 

C = 6DR B + (D B - n)SR. (3.2) 

While computing the variation of the ES equations, what we encounter most often is the 
variation (5(LJ n i£), which thanks to the ansatz (2.8) sums up to 

pin _ r n pin _ r n 

5(n n R) = ^ 1 -5UR B + U n B 8R= H + r?SR, (3.3) 

□ B - n - u B - n 

and this will turn out to be the key simplification 9 . This can be used to obtain 

S(T(D)R) = ^1EM_Z^2 C + jr lSR = (p B _ ri ) H + jr 1§R (3.4) 
□ b - ri 

where we have introduced 

Note that, expanding J~(0) in the Taylor series, one can see that 5(J-'(0)R) has no pole. 
Also thanks to J" (1) (ri) = 0, Z(D) is an analytic function. Similarly, we obtain 

S(pT(p)R) = n B (D B - ri)H + + nSR), (3.5) 
= (d"Ed v R B + d»R B d u Z) , (3.6) 

SJCx = 2E ( ri R B + r 2 ) + R B {U B - n)H - -dR^ - f ) , (3.7) 

n 

where we have used conditions Eq. (2.13), and in particular F'{r\) = 0, to obtain Eq. (3.7). 



9 Similar relationships are useful for the perturbation analysis in models with non-local scalar field as shown 
in Ref. [391. 
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3.2 Perturbing the trace equation 

Assuming <5Tm^ = 0, which is true for radiation, we get the trace equation (2.7) for pertur- 
bations linearized about a background solution: 

M 2 P bR - A<5/C£ - 2A5/Ci - 6\5{aT(a)R) = 0. (3.8) 

Substituting Eqs. (3.5) and (3.7), and using condition A\ = to eliminate Mj,, we obtain 

d^-d^Re + nR B E + 2r 2 ~ + R B U B E + 3[D B (D B - n)H + JiC] = 0. (3.9) 

This is a very important equation because it gives us information about £. It can be rewritten 
as 

V( = 0, (3.10) 

where 

V = [d»R B d, L + 2(nR B + r 2 ))Z(D B ) + W(D B ) (3.11) 

with 

W(D B ) = aF(Dfl) + (i? B + 3n)^°^^ . 

As we see, operator V depends not only on d'Alembertian but also on partial derivatives and 
has non-constant coefficients. Nevertheless, we have been able to find a non-local differential 
equation for one of the perturbative variable £, which one can use different techniques to solve. 
The form of the perturbed trace equation demonstrates the importance of the function £. 

3.3 Scalar Perturbations 

The metric perturbations can be divided into 4 scalar, 4 vector and 2 tensor degrees of 
freedom, according to their transformation properties with respect to the three spatial co- 
ordinates on the constant-time hypersurface. Different types of perturbations do not mix 
at the first order [46]- [50]. In this paper, we are going to focus on the scalar perturbations 
leaving the detailed analysis of vector and tensor modes for separate projects. Scalar metric 
perturbations are given by 4 arbitrary scalar functions 0(r, x a ), /3(r, x a ), ip(T,x a ), 7(1", x a ) 
in the following way 

ds 2 = a(r) 2 [-(1 + 2cp)dr 2 - 2d i f3dTdx i + ({1 - 2?/>}<% + 2d i d j ~ / )dx i dx j ] , (3.12) 

where r is the conformal time related to the cosmic one as a(r)dr = dt. To avoid gauge 
fixing issues we use gauge-invariant variables. There exist two independent gauge-invariant 
variables (the Bardeen potentials), which fully determine the scalar perturbations of the 
metric tensor [47]- [50] 

$ = 6 - -(a-d)' = 4> - x, * = ib + Uti = ip + Hx, (3.13) 
a 

where \ = a P + a 2 7> $ = P + Y, %(t) = a' /a, and differentiation with respect to the 
conformal time r is denoted by a prime. 

Explicit calculations show that £ is a gauge-invariant function, and it is therefore not 
possible to choose such gauge that £ = 0. Now, the function £ can be written as 



C = + 3*') - 2( n R B + r 2 )$ + (D B - n) SRqx. (3.14) 



a- 
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where SRqi is the gauge invariant part of the curvature variation: 

(3.15) 



2 



a', „o!' „ „ „ „ a! 



k 2 ($ - 2^) - 3-$' - 6—$ - 3*" - 9-*' 
a a a 



and we have used 

*□ 1 



a 2 



2$ f$ + 2^9^ + (V + 3^') 9 T 



(3.16) 



Hereafter is the comoving wavenumber. The formula (3.10) thus provides us with one 
equation for and 

It is essential to get another relationship among the Bardeen potentials and one can 
accomplish this task by considering the (i ^ j) components of the perturbed ES equations, 
Eq. (2.4). Trivially, all the terms with <JJ go away, SlCj = and <5r^ n = 0, and only few 
non-trivial pieces survive. After some algebra we get the following equation 

(n B -ri)Z(n B )( + T 1 [5R G1 + (R B + 3r 1 )($-y)} = 0. (3.17) 

This provides us with the second equation which we need in order to solve for the Bardeen 
potentials. 

It is curious to note that for radiative sources we were able to decouple the evolution 
equations for the the metric perturbations from the fluctuations in fluid. For the sake of 
completeness we have provided the equations that determine the fluid fluctuations from the 
Bardeen potentials in the appendix. 

3.4 Perturbation summary 

In this section we are able to summarize two differential equations for the two Bardeen 
potentials, <& and They read as follows: 

[dPRudu + 2(nR B + r 2 )]Z(D)C + W(D)C = (3.18) 
(D B - ri)Z(D)( + Ji [5R GI + (R B + 3n)($ - *)] = (3.19) 



where 



C = 6nR B + (p B - n)5R 



GI 

□b = — - 2^d T K = -<9 t - 3-5* k 

o z or a A a a A 

5u = ^d 2 + %$d T + \(& + W)d T 

a z a 6 a z 



a', „a!' „ T „ „ a' 



5R G1 = — k 2 {$> - 2tf) - 3-$' - 6—$ - 3^" - 9-*' 



a 2 \ a 



a' k 2 
6Db^ - 2R B $ - 6^($' + V) + 2^($ + 



a 3 a 2 



W(U) = 3T(D) + (R B + 3r x ) Jr{U) ~ ~ F 



□ - T\ 



(□-n) 2 

Clearly, equations (3.18) and (3.19) provide us with a coupled system of equations, albeit 
non-local, for the two Bardeen potentials, <3?, ^ . In our analysis, we never really need the 
quantities describing the matter fluctuations, but they can indeed be obtained by looking at 
the other components of the ES equations. We provide some details of that in the appendix. 
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4 Perturbations at late-times in the de Sitter limit 



Evidently it is a rather difficult task to solve the evolution equations for the Bardeen poten- 
tials for the most general case, even if £ vanish 10 . We can gain however significant insight into 
how the perturbations evolve within certain approximations. In particular, we are interested 
in addressing two important questions: (i) do we recover the GR limit during the inflationary 
phase, this would ensure that the inflationary mechanism of generating near scale-invariant 
perturbations can be trusted despite the higher derivative modifications to gravity, and (ii) 
we also want to see whether any pathologies develop around the bounce, this would lead to 
robustness of the bouncing mechanism, which several other non-singular alternatives seem 
to lack. In this section we will address the first question. 

4.1 Localizing equations 

At later times, t — > +oo, solution (2.19) tends to a de Sitter space-time, with R — > 4A. 
This simplifies the perturbation equations (3.10) and (3.17) remarkably. At large times 
R —7- —T2/r\, and therefore we get 

C « (□ - r^SRaj, 7>(n)«W(D)|ds 

where the subscript 'dS' indicates that we substitute Rb = — t^It\. Resulting perturbation 
equations are as follows: 

w(D)| dS (n-n)^Gi = o, (4.i) 

F(0)5R G1 + (-r 2 /n + 3ri)($ - *) = 0. (4.2) 

In this limit W is an analytic function with constant coefficients. We can therefore solve the 
first equation assuming that we know the roots of the algebraic (or transcendental) equation 
W(w 2 )|ds(w 2 — r i) = [37]. In the case of simple roots Wj this reduces to 11 

5R G i = or (□ - u;f)5R G i = (4.3) 

There is always one root uj\ = r\ but there could be other roots of W(cj 2 ), ujf with i = 2..N. 
Once 5R G i is found to be, say 6Ri corresponding to the root cof, we can express SR in terms 
of $ and and solve the equation 

SR G1 = SRi (4.4) 

along with (4.2). The most general solution is obviously an arbitrary superposition of all the 
solutions. 

Before proceeding further a couple of remarks are in order. Firstly, the fact that so- 
lutions to non-local equations can be approximated as a linear superposition of solutions to 
a collection of local equations is a rather interesting feature of these models, but it is not 
surprising and has been discussed in the literature. Secondly, the fact that there always exists 
a root = r% is also not surprising. A straight forward application of the results obtained 
in [14] tells us that for actions of the form (2.1), there is always a scalar degree of freedom, 
the famous Brans-Dicke theory belongs to this class as well, r\ is simply the mass squared 
of this scalar. Finally, we note that the roots, their number and values depend on J 7 , and in 
particular one can choose J-'s such that there are no roots other than r±. 

10 There have been steady progress in trying to solve and understand some of the properties of such infinite 
differential equations, see for instance [51]. 

n If some roots u>i are multiple roots they produce equations of the form (□ — uj^) mi 8 Rq\ — 0, where rrii is 
the multiplicity, and these must be treated differently. For the case of double roots see for instance [38]. For 
simplicity we assume that only simple roots are present. 
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4.2 Trivial case: 5R G i = 
Trivial possibility of Eq. (4.3) implies 

$ - * = (4.5) 

as it is in GR. What we then have to solve is 

/ a' k 2 \ 

5Rq = 2 [?,n B <S> - R b § + 2^<S> =0 (4.6) 
\ a 6 a z J 

Substituting the appropriate background functions for de Sitter, and Rb = 12H 2 , we can 
solve the above equation in terms of elementary functions 

<3? = r/[ci(cos(r/) + r\ sin(r/)) + C2(— f? cos(ry) + sin(7/))] , (4.7) 

where 

and we remind the readers that the inflationary scale factor in conformal time is given by 

air) = — with r G (-oo, 0) (4.9) 

Ht 

At large times clearly r/ — > as the scale factor is increasing exponentially, and all the terms 
die out. This is precisely the mode that one would obtain in ordinary GR in the presence 
of a cosmological constant and radiative sources. Within a few e- folds, the radiation density 
and it's fluctuations become irrelevant and the quantum vacuum fluctuations of the inflaton 
becomes the dominant source of metric perturbations. 

4.3 Non-trivial case, SRgi / 

We now consider solutions to the equation 

(n B -u%)5Ri = 0. (4.10) 

Acting — uj 2 on (4.2) we get 

(□ B -o; J 2 )($-^) = (4.11) 
Exact solutions to the above equation are known in terms of the conformal time: 

<S> k -* k = H(- T )l[d lk 4(-kT) + d 2k Y 2 {-kr)) (4.12) 
where J u , Y v are Bessel functions of the first and the second kind respectively, and 

,,2 . 



9 w? 



H 2 



(4.13) 



Notice that when all di k , c^fc's vanish we recover the well known GR limit $ — ^ = 0. 

We are primarily interested in the late time super-Hubble limit, k\r\ ^ 1, of these 
modes. The asymptotic forms of the Bessel functions only depend on the real part of u, lets 
call it i>r. Then we have 

lim^o Jv(x) ~ x UR and (4-14) 

i- N [ x~ VR for ur > oir\ 

lim -° r ^)~{ln(*) for^ = ( 415 ) 

We are now in a position to identify the different cases: 
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• If oj 2 < 0, then v 2 is positive, and u > 3/2. This means the Y v function grows faster 
that (— t)~ 3 / 2 as r — > 0, which in turn implies that 3> — grows at late times making 
the inflationary trajectory unstable, and the perturbations unviable for inflationary 
cosmology. 

• If u) 2 ~ 0. This is a rather interesting and special case corresponding to the presence 
of a massless degree of freedom. As it is clear from the asymptotics, in this case 
$ — ^/ tends to a constant at late times. In fact, if one assumes the usual Bunch- 
Davis type subHubble vacuum initial conditions for this field, then one will end up 
with a scale invariant spectrum for these modes. This is very similar to the curvaton 
scenario [4] that is routinely considered in inflationary cosmology and we plan to study 
the phenomenological implications of this mode in more details in the future. 

• If 9H 2 /4 > uf > 0, then v 2 is positive and < v < 3/2. Firstly, this means that the 
J v function decays at large times. The Y v function does grow but at a rate slower than 
(— t)~ 3 / 2 , so that $ — ^ decays as r — > leading to the GR relation (4.5). Inflationary 
cosmology is thus safe from these modes. 

• If 9H 2 /A < uf, then v 2 < 0, and v is imaginary, i.e. i>r = 0. The leading behavior 
of the two modes in <3? — ^ are then given by (— t)~ 3 / 2 and (— r)~ 3 / 2 ln(— r) and they 
both vanish as r — > CP . 

• If oj 2 is complex, then v 2 is also complex and whether these modes are dangerous or not 
depends on the value of vr. According to the discussion above, for phenomenological 
consistency of the model we must have vr < 3/2. 

Thus to summarize, our criteria for preserving the inflationary mechanism is that all 
the roots, Wj, should be such that the corresponding 

\u R \ < 3/2 . (4.16) 

For real roots, this simply means that they must also be non-negative. In particular since 
oj 2 = n is always a root, we must have r± > 0. As discussed in section 2.3, this is consistent 
with the observations in [13], where it was found that when r\ < 0, the de Sitter solution 
is no longer a stable solution, but it is when n > 0. Therefore, there exists a large class of 
models which are compatible with inflationary predictions. 

5 Perturbations across the bounce 

Having determined the late time behavior of the different perturbation modes around our 
inflationary bouncing solution, our next task is to see how they behave near the bounce. 
Since we have already provided the criteria to ensure that none of the fluctuations grow 
during inflation, all we now have to check is that there are no modes which become singular 
at any finite time. In particular we will focus near the bounce, as within a few e-folds of 
expansion the scale-factor is expected to resemble the de Sitter phase which we have already 
analyzed, at least this is what happens in the hyperbolic cosine bounce (2.19). In fact, for 
most part we will restrict to this special case to keep the analysis tractable. 
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5.1 Bounce Limit 

Near the bounce, even though the equations can again be simplified considerably it is a more 
complicated limit. The main difference with the previous section is that we can approximate 
Rb by a constant but cannot assume that all its derivatives vanish. We are however, going 
to look at symmetric bouncing solutions, so that near the bounce we can make the following 
assumptions about the background quantities: 

Rb ~ const = R b , Rb ~ 0, ORb ~ const = R2, 

R b (5-1) 
H = 0, H ~ const = — , a ~ 1, a = 

6 

where the subscript b stands for the bounce and numbers indicate the number of derivatives. 
Under these assumptions 

C « -2R 2 <$> + {D B - n)5R G1 (5.2) 
and equations (3.10) and (3.17) simplify to 

{2R 2 Z{U) + W{U)\ b )C = (5.3) 
(n B -ri)Z(n)( + Fi[SR G i + {R b + 3n) ($-*)] = (5.4) 

Note that even if R 2 = equations do not become identical to the ones in the previous 
section because non-local operators become slightly different due to the new constant R b 
instead of i?ds • It is more difficult do disentangle the above equations but on the other hand 
the d'Alembertian operator is simpler, thanks to the above mentioned assumptions: 

□ « -d 2 t - k 2 (5.5) 

and 5Rqi becomes a bit shorter 

SRgi = 6D b ^ -2i? 6 $ + 2k 2 ($ + 10 (5.6) 

As the first step we solve the homogeneous equation for £ (5.3). Then we express $ from (5.4) 
since the only differential operator acts on 1 there, and as the last step substitute everything 
in (5.2). The result will be an equation for 1 of the fourth order. The main question here is 
whether there are modes blowing up exactly at the bounce. 

In what follows in this Section we are going mainly to specialize to the exact solution 
(2.19) so that we can use some nice relations among the background quantities. These kinds 
of relations are not unique to this solution and the results can be easily generalized. Such 
general considerations however do not shed more light on the problem and the most important 
features are already captured in this particular solution. 

In order to consider the general bouncing solutions, and not just the hyperbolic cosine 
bounce, we must avoid the specific background relations for particular solutions. Although 
tedious, it is certainly possible to do so. Indeed, we just need to keep values R b and R 2 
general, although they are always related to each other, thanks to our general ansatz, as 

R 2 = ri R b + r 2 (5.7) 

We stress that the parameters r\ and r 2 are not initially constrained. The only reason why 
we are not elaborating on this general case is that this would significantly complicate all the 
formulae and would not provide any new insights into the problem. There is nothing special 
about the hyperbolic cosine bounce except the fact that it has a closed analytic form and 
also obeys some very elegant relations among its parameters. 
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5.2 Trivial case, ( = 

The simplest solution to Eq. (5.3) is £ = 0. This simplifies Eq. (5.4) to 

6D B ^ - 2R b $ + 2/c 2 ($ + *) + (J% + 3n)($ - *) = (5.8) 

and we can deduce $ from here as 

6D B ^ + 2/c 2 # - (R b + 3ri)* = {R b - 3n - 2/c 2 )$ (5.9) 

It is important to note that from Section 2.3, see Eq. (2.19), R b = 3n, and therefore equations 
can be separated in the small k 2 limit. Specializing to the known exact solution we have for 
k 2 = 

* + n* = o (5.io) 

Having n > we get bounded oscillations across the bounce. Then Eq. (5.3) reads 

<9 2 $ + 2ri$ = (5.11) 

where we have used R2 = — 3r 2 . These are again two bounded oscillating modes. 

If k 2 > things are more messy. Again discussing only the exact solution from section 
2.3, we have 

6^ + 4k 2 ^ + 6ri* = 2k 2 <$> (5.12) 
Using (" given by Eq. (5.2) we yield for £ = 0: 

n(* + (h 2 + n)*) + (-^ 2 - k 2 - n)(\k 2 ^ - (* + (^fc 2 + n)*)) = (5.13) 

This equation has four solutions 

* = c±±e ±I/±t , v± = -j=^-(9n + 4A; 2 ) ± ^9~r 2 + 4A; 4 (5.14) 

where signs are independent. 

We see that all the modes pass the bounce point smoothly and in the limit k — > 0, 
modes for $ and \I/ would mix together giving the above obtained solutions in k = case. 
What is reassuring is to note that the short wavelength fluctuations (k —> 00) makes the 
solutions more and more oscillatory rather than contributing to any exponential growth of 
fluctuations. This is an encouraging sign in terms of whether these kinds of theories can lead 
to a truly UV complete theory of gravity. 

One may worry that the right hand side of eq. (5.9) may become zero for a finite value 
of k for other more general backgrounds. This however, does not change the analysis of the 
equations and just shifts the point when perturbations decouple to some finite value of the 
wavenumber k. Finally, we note that these four modes should merge with the GR mode and 
the uj 2 = r\ mode that was discussed in the previous section during the de Sitter phase. 

5.3 Non-trivial case, ( ^ 

As in the previous section we first solve a simple homogeneous equation of the form 

(□ - w 2 )C « ~(d 2 + k 2 + uj 2 )( = (5.15) 
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This equation has a solution 



Ci = z ±e 



(5.16) 



and the exact behaviour depends on the value of Wj. It is clear, nevertheless, that it passes 
the bounce point smoothly. Such a non-zero £ serves as an inhomogeneous source term in 
right hand side of the two remaining equations. 
Using Eq. (5.4), we determine $ as 



where as before we have used particular relations for the exact solution Eq. (2.19). Substi- 
tuting the results in Eq. (5.2), we get 



Extra modes related to the inhomogeneous part will be either oscillatory functions of the 
frequency determined by or resonance modes if frequencies of Q and ^ coincide. We 
however stress that even the modes which look dangerous as they may grow exponentially do 
not blow up while passing the bounce point, and therefore must merge eventually with the 
late time de Sitter modes discussed in the previous section. They may continue to grow, in 
which case they must be eliminated according to our criteria (4.16), or they must be modes 
which eventually die out. 

6 Conclusions 

In this paper we studied the perturbations of a special class of string-motivated higher deriva- 
tive extension of Einstein's gravity, which is well-known to ameliorate the Big Bang cosmolog- 
ical singularity of the FLRW metric. More specifically, our action contains the Ricci scalar 
but is non-local in the sense that it contains higher derivative terms up to all orders. It 
can yield non-singular bouncing solutions similar to that of the hyperbolic cosine bounce 
which was first found in [12]. The higher derivatives lead to weakening of the gravity in the 
UV regime which facilitates the bounce for the FLRW metric in presence of non-ghost like 
radiation and the cosmological constant. 

For the first time we analyzed equations for the scalar perturbations for such non-local 
higher derivative theories of gravity in complete generality. We provided the general non- 
local equations that the two Bardeen potentials satisfy, and studied the evolution of the 
generic fluctuations in the relativistic fluid and the scalar potentials during the bounce as 
well as at late times corresponding to the de Sitter attractor solution. Although we have 
not included the dynamics and perturbations of the inflaton field, our analysis assessed the 
robustness of the inflationary scenario in terms of perturbations. We were able to arrive at 
a simple criteria such that the modes (small and large k ) pass the bounce phase smoothly, 
and consequently decay during the de Sitter phase. This means, at least at the classical 
level there are no pathologies present during the evolution. We also noted the possibility of 
realizing a curvaton scenario given the extra scalar degree of freedom present in these models 
are light. 




(5.17) 




(5.18) 
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One of the issues we haven't discussed in our paper is anisotropic deviations from 
the bouncing background. This requires looking at homogeneous (time dependent) but 
anisotropic Bianchi I type of models. (The isotropic homogeneous mode was already consid- 
ered in [13] and issues of stability were discussed.) This is an important issue to determine 
the robustness of the bouncing mechanism as bouncing/cyclic models are notoriously plagued 
with Mixmaster type chaotic behavior related to the growth of anisotropies in a contracting 
universe. The generalization of the application of the ansatz (2.8) that has been performed 
in this paper to arbitrary metrics (potentially inhomogeneous and anisotropic) should help in 
addressing this issue, but we leave this for future considerations. For stability analysis related 
to anisotropies in other non-local gravity and string-inspired models one can see [25, 52] and 
references therein. It is probably worth pointing out though that anisotropies are expected 
to be the largest at the bounce point, and thus if there is a "smooth-enough" bouncing patch, 
anisotropies are only going to decrease at large times in the past or future. 

For simplicity, some of our analysis was done for the hyperbolic cosine bounce, but 
we discussed how the results should generalize to other bouncing solutions. This strongly 
suggests that higher derivative actions of gravity considered in this paper do indeed yield a 
non-singular bounce without a pathology, and in doing so can provide a model of inflation, 
both in the past and the future. 
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A Equation for the matter density perturbations 

Scalar stress-energy tensor perturbations are parameterized as follows 



where p is the energy density, p the pressure, 5p and 5p are there respective variations, v s the 
velocity or the flux related variable and n s the anisotropic stress. The perturbation functions 
in Ty are as follows: Sp(rj, x) = 5p(r], k)e lk i xl and similar for dp, v s and ir s . The corresponding 
often used gauge invariant quantity is [48]- [50]: 



o — — 




.0 
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In this paper, we consider only radiation for which 



TT S = 0. 



P=\p; 



-Sp. 



From the perturbed (0, 0) and (0, i) ES equations we express v s and 5p via the Bardeen 
potentials: 



5p=X 



+ (nR B + r 2 )H + 2D B (D B - n)H + 2Ji(C + n5R G i) 



-AF^Rb + 3n) (3^(*' - + ^J*) - 6^T + i? B T + ^ (T" - HT') 



and 



3A;A 
2pa? 



2F 1 {R B + 3n) - H$) + \b! b B + T' - Kt 



where 



T = (D B - n)H + Jiffloi and H = Z(D)(. 
The resulting expression for e then is 



A 



— f-^' + ( ri R B + r 2 )H + 2D B (D B - n)H + 3^H + 2Ji(C + nfflci) 



4-2^i (i? B + 3n)* + — ^.R^S + <^ i? B - 12-2- ^ T + -=T" 



(A.l) 



(A.2) 



(A.3) 



One interesting special case is R B — > —Sri, but this is not a limit for the known background 
solutions. 
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